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Abstract. Pinpointing the progenitors of long duration gamma-ray bursts (LGRBs) remains an 
extremely important question, although it is now clear that at least a fraction of LGRBs originate in 
the core collapse of massive stars in type Ic supernovae, the pathways to the production of these 
stars, and their initial masses, remain uncertain. Rotation is thought to be vital in the creation 
C"| | of LGRBs, and it is likely that black hole creation is also necessary. We suggest that these two 

O |. constraints can be met if the GRB progenitors are very massive stars > 20 M© and are formed 

in tight binary systems. Using simple models we compare the predictions of this scenario with 
observations and find that the location of GRBs on their host galaxies are suggestive of main- 
sequence masses in excess of 20 M©, while 50% of the known compact binary systems may have 
been sufficiently close to have had the necessary rotation rates for GRB creation. Thus, massive 
i— i , stars in compact binaries are a likely channel for at least some fraction of LGRBs. 
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INTRODUCTION 

The link between some long-duration gamma-ray bursts (LGRBs) and stellar collapse 
is now firmly established [1, 2, 3] In particular LGRBs appear to originate in type Ic 
supernovae and frequently in hypernovae. There is significant diversity in the LGRB 
population, including the highly luminous ("classical") burst population, originating 
from a mean redshift of 2.8 [4], a nearby lower luminosity sub-class, which has a much 
higher space density (e.g. [5]), and possibly some LGRBs which are not associated with 
supernovae. This suggests there may be some variety in the progenitors of the bursts, 
and that they may not come from systems undergoing identical evolution. Nonetheless, 
there remains a strong drive to identify the progenitor stars themselves. 

In the standard picture, core-collapse supernovae lead to LGRBs when the stellar core 
has a rapid rotation just prior to collapse. Some of the inf ailing material then forms a 
torus around the newly-formed central compact object (often assumed to be a black 
hole) and subsequent accretion of the material in the torus then fuels the LGRB. We are 
naturally left with two questions: 1) What leads to the high rotation rates required to 
form a torus? 2) What is the minimum progenitor mass required (or, in other words is 
the formation of a black hole (BH) rather than a neutron star (NS) necessary)? 

For massive single stars, it is not clear whether sufficiently high central rotation 
rates may be maintained to produce a torus on core collapse (e.g. [6, 7]), although it 



has been suggested that rapidly rotating metal-poor stars can retain sufficient angular 
momentum [8, 9] Alternatively, binary scenarios suggest a way of removing the envelope 
and providing a source of angular momentum (e.g. [10, 11]). We have explored this idea 
further [12] and present a summary of the key results here. 

In addition, we have explored whether the observation of Fruchter et al. [13] that 
LGRBs are more concentrated on their host galaxy light than SNe may be due to the 
greater mass of a LGRB progenitor than a typical SNe progenitor [14]. A review of the 
key results of this work is also presented here. 

NEUTRON STAR BINARIES AND LONG-DURATION 
GAMMA-RAY BURSTS 

Evolutionary pathways to compact object binaries 

Binaries containing two compact objects (ie white dwarfs, neutron stars or black 
holes) can be formed through several channels (see e.g. [15]). The basic scheme is 
shown in Fig. 1, and is essentially the same as that described in Bhattacharya and van 
den Heuvel [16]. 

The primary evolves first to produce either a black hole or neutron star, possibly 
following a phase of mass transfer to the secondary. When the secondary evolves, and 
fills its Roche lobe, the formation of a common envelope phase is likely owing to the 
mass ratio of the system. The black hole/neutron star and the helium-star core of the 
secondary will then spiral towards each other as the enshrouding gaseos envelope is 
ejected. The post-common envelope system will be a compact binary. A key point is that 
the helium-star core will be rotating rapidly because of tidal locking. In at least some 
systems, the rotation rate will be sufficiently large that when the secondary explodes as 
a core-collapse supernova, some of the infalling material will form a torus around the 
newly-formed neutron star or black hole. 



Discs around neutron stars and black holes 

In a core-collapse supernova, the material within the stellar core collapses to form 
a black hole or neutron star. If a core of a massive star is rotating sufficiently rapidly, 
some of the infalling material will be centrifugally supported (as it will carry its angular 
momentum with it) and form a torus around the central black hole or neutron star. The 
torus will have a radius of DGM c /c 2 (where M c is the mass of the black hole or neutron 
star) if the specific angular momentum of the infalling material is given by \/DGM c /c. 
In order to have a stable orbit around a (non-rotating) black hole, we require D > 6 (or a 
radius of around 12 km for a 1.4 M Q black hole). Forming a disc around a proto-neutron 
star requires more angular momentum as a disc would have to be formed with a radius 
somewhat larger than 10 km. 

If we assume that tidal locking occurs at the beginning of the helium main sequence 
and that for the late stages of evolution of the helium star the core decouples from the 
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FIGURE 1. An evolutionary pathway to the creation of a binary containing a rapidly-rotating core- 
collapse supernova in a tight orbit. The primary evolves first, possibly transferring material to the sec- 
ondary (stage 1). It then produces a neutron star (NS) or black hole (BH), when it explodes as a core- 
collapse supernova (stage 2). The secondary then evolves, filling its Roche lobe (stage 3) and transferring 
material to the NS/BH producing a common envelope phase (stage 4). The NS/BH and He core of the 
secondary spiral together ejecting the surrounding envelope producing a very compact binary (stage 5). 
Tidal locking produces a rapidly-rotating He star such that the rotation is significant when the secondary 
explodes as a core-collapse supernova, with a torus being formed around the central compact object by 
infalling material. 



envelope. In this case, following Podsiadlowski et al. [11] we can compare the required 
angular momentum at the time of collapse with that at the edge of the iron core at the 
start of the helium main sequence, and, assuming tidal locking equate this to an orbital 
frequency (i.e. (0 = \fDGM c /R^.c, where R c is the radius of the iron core). Assuming 
synchronous rotation, this gives a critical orbital separation [12] 

a < (4M tot c 2 ^/9DGM c 2 ) 1/3 , (1) 



a\ 60 (R C \ A/3 (M ioi \ l " (M^- 2 ' 3 



or 

where M c and R c are the core mass and radius, and M tot is the total mass of the binary. 
We use here the models of Heger et al. [17], the key parameters are M c = 1.7 M , 
R c ~ OARq, for a He star of mass 7.71 M , this yields a critical separation of ~ 3 R 
for the formation of a disc at the innermost stable orbit of a 1.7 M Q BH. 

Do systems exist which are sufficiently tight to form discs when the secondary un- 
dergoes a core-collapse supernova? In order to answer this question, one can consider 
the observed compact binaries containing black holes or neutron stars which would be 
formed after the second supernova explosion (assuming the binary remained bound). 
There are currently no known black hole-neutron star binaries. However there are cur- 
rently eight NS-NS binaries and three WD-NS binaries. The semi-major axis and ec- 
centricity of a binary provides limits to the separation of the binary at the moment of 
the second supernova. In Fig. 2 we plot the semi-major axes and eccentricities of the 
observed systems. The plot also illustrates the minimum separation required to form a 
disc around the second compact object for discs of various radii. From this figure, we 
can see that discs are likely to have formed during the second supernova in about half 
of all osberved binaries. Energy released from such discs may power some of the long 
gamma-ray bursts, as discussed below. 



The contribution of compact binaries to long gamma-ray bursts 

Long-duration gamma-ray bursts are commonly thought to originate from discs 
around black holes located at the centre of a supernova. However, one can also con- 
sider whether discs around newly-formed neutron stars may also produce some form of 
gamma-ray burst. The maximum energy released in the accretion of the torus is given 
by E = GM ns/bh M acc /fl ns/bh , or £ ns = 3.6 x lO 53 (M acc /M ) ergs, for a 1.4 M NS and 

£bh = 1 x lO 54 (M acc /M ) ergs, for any mass BH (since Rbh)- 

The extrapolation from disc accretion to gamma-ray luminosity is far from trivial 
since it requires an assumption about the conversion of accretion luminosity into y-ray 
energy. One plausible mechanism of providing this energy is via neutrino-antineutrino 
annihilation. If this is assumed to be the energy source then the accretion energy can be 
related to the observed gamma-ray energy via several efficiency factors which account 
for the conversion of accretion energy to neutrinos, the cross section for neutrino - 
antineutrino annihilation, the subsequent fraction of energy that is transferred into a 
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FIGURE 2. The eccentricity and semi-major axis of observed systems containing either two neutron 
stars, or a neutron star and white dwarf. The error bars give an indication of the range of separations 
between the two stars during their orbit. Neglecting the inspiral the separation at the time of the supernova 
must be taken from this range (see Table 1 for the calculations including the effects of inspiral via 
gravitational radiation). The three vertical lines represent the critical separations necessary at the time 
of core collapse for a centrifugally supported disc to have formed at a distance of 6M, 20M and 50M from 
the newly formed compact object, based on the models of Heger et al. [17] described in section 3 (where 
M = GM„ s / hh /c 2 and is — 12,40 and 100 km for D=6,20 and 50 for a 1.4 M NS). Lower total masses, 
as were likely the case of J0737-3039 require slightly tighter orbits, although the orbit is only a weak 

1 /3 

function of the total mass (M U)t ). If the binary separation is less than this (i.e. to the left of the line) then 
disc formation is favoured. The data are from Champion et al. [18] and Lorimer [19]. In cases where the 
masses of the two components have not be measured 1.4 M has been assumed. Known NS-NS binaries 
are indicated with circles, while NS-WD or those with uncertain companions are marked with open circles 



baryon free jet, and finally the fraction of this energy which is emitted as gamma- 
rays [20]. Following Oechslin and Janka [20] we assume that the product of the these 
efficiencies is ~ 10 3 . Thus the observed luminosities of low luminosity GRBs of 
10 48 - 10 50 ergs can be explained by the accretion of 0.01 < (M acc /M Q ) < 0.3 of 
material from the disc. 

More-massive versions of the systems which formed the observed NS-NS binaries are 
expected to form BH-BH binaries and are thus candidates for classical long gamma-ray 
bursts. Could both supernovae in such systems produce a gamma-ray burst? One can 
consider the separation required in order to give the core of the primary a sufficiently- 



high rotation rate. It is found that the size of the stellar is restricted to less than a few solar 
radii. In other words, the primary of a binary containing two massive stars is unlikely to 
be tidally spun-up sufficiently to produce a gamma-ray burst. However the black-hole it 
produces will work effectively to tidally spin-up the core of the secondary star following 
a pathway similar to that shown in Fig. 1 . 

A NEW CONSTRAINT FOR GAMMA-RAY BURST PROGENITOR 

MASS 

The basic idea 

Recently Fruchter et al. [13] have conducted a survey of the galactic environments of 
both long duration GRBs and core collapse SNe (i.e. all types of core collapse events, 
including SN II, lb and Ic). These results demonstrate that GRBs are highly concentrated 
on their host light, significantly more so than the core collapse supernova population. 
Fruchter et al. [13] further suggest that this can be explained as being due to the 
GRBs originating in the most massive stars, which, upon core collapse form black holes 
rather than neutron stars. Here we further explore this possibility and attempt to derive 
plausible limits on the progenitor lifetime and mass based on the observed distributions 
of core collapse SNe and GRBs upon their host light. Using a simple model, motivated 
by the distributions of young star clusters in a local starburst galaxy, we explore the 
expected distributions of stars of different masses upon their host galaxies and compare 
these to the observed distributions from Fruchter et al. [13]. Our results demonstrate 
that for plausible models more massive stars are always more concentrated on their host 
light than lower mass stars. Further, given that supernovae originate from stars with 
initial masses > 8 M , we find that the observed distributions of GRBs on their host 
galaxies can naturally be explained by progenitors with initial masses in excess of 20 
M . 



The model 

GRB host galaxies at high redshift are typically starburst galaxies. The Antennae 
(NGC 4038/4039) are a natural, local, analogue and were used as a template for con- 
structing a simple model for GRB host galaxies. In this model we consider the popula- 
tion of young stellar clusters, as well as the stellar population contributing to the back- 
ground light within the galaxy. The galaxy model consists of several key parameters: the 
surface density of clusters; the distribution of cluster masses; the distribution of cluster 
ages; the distribution of background light; and the distribution of clusters on the back- 
ground light. The first two properties were taken from observations of NGC 4038/4039 
when viewed at z ~ 1. The surface density of clusters expressed in terms of number 
of clusters per pixel is ~ 0.15, although this is far from uniform across the galaxy. We 
use the observed distribution of cluster masses, which follows dN /dM c \ M d 2 , with 
^ci,min = 4 • 10 4 M Q and M cl max = 10 6 M Q , where M cl is the cluster mass. The age dis- 
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FIGURE 3. Fraction of objects plotted against fraction of light for observed SNe (red) and GRBs (blue) 
together with the results from our model (black lines). Black lines from top to bottom correspond to 
minimum progenitor masses of 8, 20, 40, 60, and 80 M . 

tribution of clusters was taken from Fall et al. [21] and follows dN /dx T -1 , where 
1 is the cluster age. In order to mimic an (almost) instantaneous burst of star forma- 
tion, clusters are created in the model according to this distribution over a period of 10 7 
years. The distribution of the background light adopted lies in the middle of the distri- 
butions observed for galaxies hosting supernovae and gamma-ray bursts and is given by 
dN/dLpa L pi J 5 with L piXjmax /L pixMn = 20, where L pix is the luminosity of a pixel. 
We distribute the stellar clusters such that there is a correlation between background 
brightness and the number of clusters present in a pixel. The degree of correlation was 
chosen to give a good match to the cluster distribution observed in NGC 4038/39. 



Results 

Using the parameters described in the previous section we performed runs for min- 
imum progenitor masses of 8, 20, 40, 60, and 80 M . The results are shown as black 
lines in Fig. 3 together with the observed distributions of SNe (in red) and GRBs (in 
blue) from Fruchter et al. [13]. 
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FIGURE 4. KS-probabilities of our model results following the observed SN (red) and GRB (blue) 
distributions. The probabilities are plotted as a function of the minimum progenitor mass and were 
calculated for masses of 8, 20, 40, 60, and 80 M . A spline has been fitted through the data points. 



The model distributions for all masses were KS-tested against the observed SN and 
GRB distributions and the resulting probabilities are shown as a function of mass in Fig. 
4. While the probability of following the SN distribution decreases with increasing mass, 
the likelihood of following the observed GRB distribution increases rapidly from 8 to 40 
Mq and then flattens out, reaching a weak maximum around 60 M & . The shapes of the 
two probability functions look the same for all realisations of the model, although the 
peak probabilities can change by about 0. 1 between different runs. These results strongly 
suggest that GRB progenitors are significantly more massive than SN progenitors. 

SUMMARY 

We have suggested that GRBs originate from the most massive stars [14] and that many 
of these stars exist within tight binary systems [12]. After the GRB a NS-NS, NS-BH 
or BH-BH binary may remain, as the separation at the time of the GRB is moderately 
small many of the systems will merge within a Hubble time and might be observable as 
a short GRB. These systems are spectacular, at some point in their evolution they will 
be observed as two supernova, two GRBs as well as X-ray binary (and possibly even 



ultraluminous X-ray source) phases. 

We note that while it is likely that LGRBs are formed via this channel it is not 
necessarily the only one operating, indeed the diversity in the properties of observed 
GRBs implies there may be significant variations in the progenitors themselves, for 
example the low luminosity burst population may be formed via neutron star creating 
supernovae [12, 22] while single stars, which undergo complete mixing on the main 
sequence (e.g. [9]) may also create a subset of the GRBs. 
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